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^ I Monte Carlo simulations of the Asakura-Oosawa (AO) model for colloid-polymer 

'^ ■ mixtures confined between two parallel repulsive structureless walls are presented 

O . and analyzed in the light of current theories on capillary condensation and interface 

I localization transitions. Choosing a polymer to colloid size ratio of g = 0.8 and 

^ ' studying ultrathin films in the range of Z) = 3 to Z? = 10 colloid diameters thick- 

ly ! 
0^ . ness, grand canonical Monte Carlo methods are used; phase transitions are analyzed 

'^ i via finite size scaling, as in previous work on bulk systems and under confinement 

O ' between identical types of walls. Unlike the latter work, inequivalent walls are used 

oo : 

^D ' here: while the left wall has a hard-core repulsion for both polymers and colloids, 



X 



at the right wall an additional square-well repulsion of variable strength acting only 



H ' on the colloids is present. We study how the phase separation into colloid-rich and 

colloid-poor phases occurring already in the bulk is modified by such a confine- 
ment. When the asymmetry of the wall-colloid interaction increases, the character 
of the transition smoothly changes from capillary condensation-type to interface 
localization-type. The critical behavior of these transitions is discussed, as well as 
the colloid and polymer density profiles across the film in the various phases, and the 
correlation of interfacial fluctuations in the direction parallel to the confining walls. 
The experimental observability of these phenomena also is briefly discussed. 



I. INTRODUCTION AND OVERVIEW 



When fluid systems are confined in nanoscopic pores or channels, one expects that the 
pha^e behavo. »„ be profoundly .odmed J fl fl fl fl fl B. B fl Q. Sucb effec. 
have found an increasing attention recently, for instance because of the current interest 
to fabricate devices of nanoscopic size and to manipulate chemical reactions in nanoscopic 

jlll . Il2l . I13I . 1 14 . I15I . Ilq . In addition, porous materials 



reaction volumes ("lab on a chip"), etc. 



with pores of nanoscopic widths are useful as catalysts or for applications such as mixture 
separation, pollution control, etc. p, [iT], [isl, [l9[ . 

However, such applications often are based on empirical knowledge, the theoretical un- 

10|. In order to make 



derstanding of confined fluids still being rather limited p, [Zl, M, M-, 
progress with the theoretical description of fluids under confinement by the methods of 
statistical thermodynamics, it is desirable to start with relatively simple model systems, 
where both the geometry of confinement is well characterized, and the relevant interactions 
among the fluid particles and between the fluid particles and the confining solid surfaces 
are sufficiently well understood. Last but not least, suitable experimental tools should be in 
principle available to put the theoretical predictions to a test 
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For these purposes it is hence useful to consider colloidal suspensions [20 
24l |. exploiting the analogy between colloidal fluids and fluids formed from small molecules, 
but taking advantage of the much larger length scales (in the /im range), of the colloidal 
particles. Such systems allow detailed experiments in which individual particles can be 
tracked through space in real time using confocal microscopy techniques 25|]. Particularly 
useful systems in the present context are colloid-polymer mixtures, which can undergo in 
the bulk a liquid- vapor like phase separation into a colloid-rich phase (the "liquid" ) and a 
colloid-poor phase (the "vapor") 23|, |26|. This phase separation is due to the (entropic) 
depletion attraction between the colloids caused by the polymers. A very simple model. 



the bulk 
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due to Asakura and Oosawa [27|] and Vrij [28|] describes the resulting phase separation in 



331] in excellent qualitative agreement with the experiment 23 1. 



While initially it was thought that mean fie 
Monte Carlo (MC) simulation resu 
extensive MC simulation study 32, 



ts 



30 



d theory 29|] accounts very accurately for the 



311 ] of this Asakura-Oosawa (AO) model, a more 



331 ] revealed clear evidence for Ising-like critical behavior 



341 ] over a broad regime of control parameters. 



When such a colloid-polymer mixture is confined by hard walls, also a depletion attraction 



of the colloids and the walls occurs 35|] and can cause (in semi-infinite geometry 36 
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401 ) the formation of wetting layers 



interfacial tension between unmixed phases 47, 



wave fiuctuations 
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37|, 



461]. Due to the very low 



50(1, thermally activated capillary 



55l | are readily observable in experiment 56| and simulation 



[50|. The phase behavior of colloid-polymer mixtures in confinement can be also studied 
experimentally. Therefore, this issue has been addressed in recent computer simulation 
studies, considering the confinement of colloid-polymer mixtures by two parallel hard walls 



a distance D apart 



um 
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59| . These studies have confirmed the fact that lateral phase 



separation in a thin film geometry exhibits a critical behavior belonging to the class of the 



two-dimensional Ising model 58|. Also the the scaling relations of Fisher and Nakanishi 60] 
have been verified. Unlike the case of confinement of small molecule fiuids in nanopores, the 
size of the particles in colloidal fiuids by far exceeds the scale of the atomistic corrugation 
of the pore walls, and hence the effects of this corrugation on the packing of particles near 



the walls 



61 



621] need not be considered here. 



A very useful aspect of colloidal suspensions is that interactions among such particles can 



be tuned by suitable surface treatment 20|, 
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631 ]. E.g., a short-range repulsion between 



641 ]. Similarly, 



colloidal particles often is created by coating them with a polymer brush 
one could cancel (partially or completely) the depletion attraction of colloids towards a hard 
wall by coating the latter with a polymer brush, choosing the grafting density and chain 
length of these fiexible polymers appropriately. In a colloid-polymer mixture, however, 
for moderate chain stretching in the polymer brush the polymers in the solution still can 
penetrate into the brush, experiencing hence a much weaker interaction than the colloidal 
)articles. Only for strongly stretched chains, as occurring in very dense polymer brushes 



651 ]. a repulsion of the polymer coils in the solution would result as well, even if the chemical 



nature of the polymers in the solution and in the brush is identical ( "autophobicity effect" 



67|) 



This tunability of the wall-colloid interactions opens the possibility to realize a situation 
of a slit pore with asymmetric walls: suppose the left wall is simply a hard wall, attractive for 
the colloids, and the right wall a coated hard wall, repulsive for the colloids (Fig.[T]) 68(]. With 
a colloid-polymer mixture confined between such asymmetric walls, the possib ility arises to 



realize the "interface localization transition" 



i7i,y,y, 



Id, 
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78|. This 



transition is illustrated in Fig. [H Here, the so-called "polymer reservoir packing fraction" 
is defined by t]^ = (47r/3)i?pexp(/ip/fcBT) (with _Rp and yUp the radius and the chemical 
potential of the polymers, respectively) and plays the role of inverse temperature when 
we compare the behavior to that of a fiuid of small molecules that undergoes a liquid- 
vapor transition. While in the bulk colloid-polymer mixture phase separation sets in when 
the variable t]^ exceeds the critical value fj^cnt^ ^^is transition is rounded in the thin film. 
Starting out from a layer enriched with colloids on the left wall and enriched with polymers 
at the right wall, a stratified domain structure forms, with a domain wall separating the 
colloid-rich phase in the left part and the polymer-rich phase in the right part of the slit 
pore (state BI in Fig. [1]). Only at a much larger value ^pcrit(-^) ^ sharp phase transition 
occurs in the thin film, with the colloid-polymer interface being bound either to the right 
wall (phase B Ila) or to the left wall (phase B lib). Along the fine fi = /icoex(-D,?7p) these 
two phases may coexist. 

Of course, in an experiment one does not have at one's disposal the intensive variables n 
and the "polymer reservoir packing fraction" rj^, but rather the volume fractions of colloids 

and polymers, 

4:71 An 

r^^ = —RlNjV, r^^ = —RlNjV, (1) 

where V is the volume of the system. Re the radius of the spherical colloidal particles, and 
Nc, Np are the particle numbers of colloids and polymers, respectively. Since rjc, rj^ are 
densities of extensive thermodynamic variables, the first order transition lines /icoex(-D, r/p 
in the plane of variables rjc, rjp are split into two phase coexistence regions. Bringing the 
thin film from the one-phase region to inside the two-phase region (e.g. by adding polymers 
to the solution), one creates a state of the slit pore where in parts of the system the interface 
is bound to the left wall and in other parts it is bound to the right wall. These phases are 
then separated by interfaces running across the film from the left to the right wall (or vice 
versa). A similar phase coexistence between the two phases AI, All occurs in the case of 
capillary condensation- like transitions for symmetric walls (left part of Fig. [1]). As always, 
the amounts of the coexisting phases is controlled by the lever rule. 

In the limit D -^ oo oi the film thickness, we recover a semi-infinite system and then 
wetting transitions are expected to occur, so that, in the symmetric wall case, in the region 



'7pcrit < Vl < Vpw ^OT fi = /icoex(oo) both walls are 



the walls are nonwet ("incomplete wetting" [Sy, 



m 



cornpletely) wet, while for rj^ > rj^ 



38, 



39 



p,w 



40]). In fact, the colloid- rich 



surface enrichment layers indicated for the phase All are the precursors of wetting layers 



that appear when D ^ oo. Of course, no (infinitely thick [Sy, l37|, l38|, |39|, |40|]) true wetting 
layer fits into a thin film of finite thickness D, and thus the wetting transition at 1]^ = rjl^^ 
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m. 



J,Q)i 



is rounded off in the thin 



(which we have assumed to be of second order 
film. 

For asymmetric walls in the limit D -^ oo the wetting transitions at both walls will occur, 
in general, for different values of 7]^ at both walls. In Fig. [1] we have arbitrarily assumed 
"f^. In the simplistic Ising model with "competing surface magnetic fields" 



;hat_^;^f^> ?]Jfjf^*. In the 

d, y , y, y, y , m hi 



and Hd, one can consider a situation with Hd = —Hi, where 
these transitions then coincide, rj^^^^ = "/^p"^^*- However, such a special symmetry never is 
expected for a colloid-polymer mixture (which has an asymmetric phase diagram already 
in the bulk). Note, however, that ioi D —>■ oo one does not expect that for interface 
localization transitions ^pcrit(-^) converges to the bulk critical point, ''7pcrit- rather one 
expects a convergence towards the wetting transition which is closest to the bulk transition 

0. 

In the present paper, we shall present evidence from Monte Carlo simulations that the 
scenario sketched in Fig. [1] is correct, and we shall characterize the behavior of colloid- 
polymer mixtures confined by asymmetric walls in detail, considerably extending preliminary 
work 68|. Extensive results for the case of symmetric walls have been presented earlier 58|, 
59|]. As in previous studies in the bulk [321, l33| the simulations are carried out mostly in the 
grand-canonical ensemble, using a dedicated grand-canonical cluster alg orithm 321] together 
with re-weighting schemes such as successive umbrella sampling 79]. Phase transitions 



are analyzed by finite size scaling methods 80, 



81 



82| . varying suitably the lateral linear 



dimensions L along the walls. For a description of these techniques, the reader should consult 



our earlier work 



58 



59|. 



In Sec. II we now present a study of the "soft mode" phase [72] BI for a relatively thick 



film (thickness D = 10 colloid diameters). Such phases with de^ 



great interest due to their large interfacial fluctuations 72 



ocalized interfaces are of 
73l . 174 l83l . l84| , and consequences 
of such fluctuations have been seen in experiments both on polymer blends 85| and colloid- 
polymer mixtures 46|]. Sec. Ill then gives a discussion of the interface localization transition 
for an ultrathin film [D = 3), attempting to verify the above statement that the critical 
exponents should be those of the two-dimensional Ising model. Sec. IV discusses the phase 
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behavior when both film thickness and the strength of the short range colloid-wall repulsion 
are varied. Finally, Sec. V summarizes some conclusions. 

II. FORMATION AND PROPERTIES OF THE INTERFACE IN THE SOFT 

MODE PHASE 

All our Monte Carlo simulations refer to the standard Asakura-Oosawa (AO) model and 
use the same size ratio q = Rp/Rc = 0.8 as the previous work in the bulk [32, l33| and for 
symmetric walls 58|, |59|. In this case, it is known that the critical point in the bulk occurs 



at 



32 



33| 



% 



p,crit 



0.766 ± 0.002, r^ccrit = 0.1340 ± 0.0002, r/p^^rit = 0.3562 ± 0.0006, 



(2) 



and also the coexistence curve between the colloid-rich phase (?7c/) and the polymer-rich 
phase {r]c,v) is known rather precisely, as well as the interfacial tension 32, l33|, |50(]. We 



now consider a L x L x D geometry, where all lengths are measured in units of the colloid 
diameter 2Rc, and periodic boundary conditions are applied in x and ?/-directions only. For 
the thickness D, the values D = 3, 5, 7, and 10 are used, while the linear dimension L in 
parallel direction is chosen in the range from L = 15 to L = 30. The left wall, located 
at 2 = 0, is taken purely repulsive for both colloids and polymers. As for the interaction 
between the colloidal particles (which is infinite if two colloids overlap and zero else, as 
well as the colloid-polymer interaction which also is infinite if a colloid particle overlaps a 
polymer and zero else), we take a hard wall repulsion. 



KAz) 



oo, z< Re, Ul^z) 



0, Z > Rr 



Ui^z) =oo, z< Rp, Ui^z) =0, z> Rp, 



(3) 
(4) 



for both colloids [U^^z)] and polymers [f^^.pl-^)]- ^^ ^^e right wall, however, we add a 
square well potential of strength e and with an additional range Re- Thus, the potential 
acting on the colloids is 



KA^) 



= 0, z < D- 2Re, 
UlAz) = ^, D-2R,<z<D-R, 
UlA^)=oo, z>D-R,. 



(5a) 
(5b) 
(5c) 



This square well potential [Eq. (15bp ] could be realized by a polymer brush of low grafting 
density and height Re, for instance, so that the region of z where the colloid penetrates into 
the brush leads to a finite energy penalty e only (note that we use the convention that the 
temperature fc^T = 1; of course, one could also consider square well potentials of arbitrary 
range). For the polymers, on the other hand, the interaction is taken to be of the same type 
as in Eq. (JH), 

Ul^Jz) = 0, z<D-Rp, Ul^^{z) = 00, z>D-Rp. (6) 

This potential models the interactions of polymers with a hard wall coated with polymer 



brushes: Under good solvent or Theta solvent conditions [86|, polymers can overlap with 
weakly stretched polymer brushes with little free energy cost. 

It turns out that a phase behavior as sketched in the right part of Fig. [TJoccurs ii e > 2.5. 
Figure [2] presents some typical profiles of the average local volume fraction of colloids ric{z) 
and polymers rip{z) across the slit pore, for the case e = 2.5 and D = 10. Panel (a) shows the 
profiles for 7]^ = 0.18 and rj^ = 0.7, corresponding to a state point where the bulk colloid- 
polymer mixture is still in the one-phase region. Nevertheless, the profiles of ric{z) and 
rip{z) exhibit pronounced inhomogeneities: the polymer profile rjclz) displays a pronounced 
peak close to the right wall, and decays with increasing distance from the right wall to a 
plateau, almost independent of z, in the regime 3 < z < 6. Very close to the left wall, 
where the volume fraction of colloids is strongly enhanced, the concentration of polymers 
is also inhomogeneous (indirectly induced by the colloids, since polymers and colloids must 
not overlap), before rjp^z) abruptly decreases to zero for z = Rp. The colloid profile ric{z) 
shows a very pronounced peak close to 2; = i?c, on the other hand, which can be attributed 
to the depletion attraction of the colloids to the hard wall. One can recognize a second peak 
near z = 1.6 and a weak third peak near z = 2.5, these peaks represent the well-known 
"layering" of hard particles near smooth repulsive walls. In the central part of the thin film, 
for 3 < z < 6, the profile ric{z) is almost flat; thus the surface enrichment of the colloidal 
particles at the hard wall is a short range effect. In the regime near the right walls, where the 
polymers are attracted, we recognize first a smooth decrease of ric{z) in the range where the 
pronounced increase of %(-z) sets in. For z = D — 2Rc = 9, where the additional repulsive 
potential sets in, a downward step in ric{z) occurs, as expected. 

It is interesting to contrast the behavior in panel (a), showing surface enrichment of 
colloids (left) and polymers (right) at the walls confining an otherwise homogeneous mixture. 
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with the behavior in panel (b), which refers to a state where in the bulk phase separation 
has occurred. Indeed, Fig. [2]d gives rather clear evidence for a phase separation in the z- 
direction perpendicular to the confining walls, of the type denoted as BI in Fig. [H The 
polymer rich phase occurs on the right side of the thin film, and ?7p(-2) reaches very small 
values for z < 4. Near z = 6 we recognize inflection points in both profiles ?7p(2;), ^/cl^) as are 
typical for interfaces between coexisting phases. Again the profile riJ^z) exhibits the typical 
layering oscillations for small z. No such layering occurs for the polymers near z = D, oi 
course, since the polymer-polymer interaction is zero, the polymer-rich phase is like a dense 
ideal gas. 

Panels 2c and 2d illustrate states corresponding to the phases Bllb and Blla in Fig. [H 
respectively. In the polymer-rich phase the interface position is at about z = 2.5, and unlike 
Fig. [2^ (where the interface is freely fluctuating in the center of the slit pore) the width of 
the interface is only about two colloid diameters. Such a state is typical for a colloid-polymer 
interface tightly bound to the left wall. Figure [2li is the counterpart showing the profiles 
in the colloid-rich phase, where almost all polymers are expelled, apart from the immediate 
neighborhood of the right wall. 

We conclude that these profiles do give qualitative evidence for the existence of all three 
phases BI, Blla and Bllb in Fig. [1] We now study the phase with the delocalized interface 
(BI) more closely. In particular, we are interested in how the interfacial profiles change when 
the inverse-temperature-like variable rji^ is varied (Fig. [3]). Defining an order parameter m 
and the coexistence diameter 6 as follows, 

m = {r^i-r^:)/2, 5 = (r/,^ + r/f)/2, (7) 

we choose the average volume fraction of the colloids such that rj^ = S, and we attempt to 
fit the colloid density profile by a tanh function, 

ric{z) = S — mtanh[{z — Zo)/w] ; (8) 

here zq is the position of the interface center and w is the interfacial width. Fig. [3^ shows 
that Eq. ([8]) provides a good fit of the colloid density profile, for all values of 1]^ from 0.90 
to 1.10. For rjp = 0.80, however, the profile is extremely wide, due to the proximity of the 
critical point in the bulk [Eq. ([T])], and then the fit is less convincing. Indeed, the polymer 
density profile rip{z), Fig. [3]d, for rj^ = 0.8 does not even exhibit an inflection point, while 



9 



for all larger values of rj^ an inflection point clearly is present (it occurs roughly aX z = zq, 
the inflection point of the polymer density profile, which is roughly at zq ^ 0.20 ± 0.05). 

Of course, one notes that ric{z) does not reach the regime of homogeneous "liquid" density 
rj^, since for z < —1.5 in Fig. [3|l the layering effect caused by the repulsive wall aX z = —5 
already sets in. Likewise, the surface enrichment of the polymers at the right wall distorts 
the profiles for z > 3.5 in Fig. [3]d. We also note that the profiles seem to have common 
intersection points (which do not coincide with zq, since both m and 6 depend on rjV). The 
common intersection point of the colloid profiles is at 2; = 0.5 ± 0.05, while the common 
intersection point of the polymer profile is at 2; = 0.0 ±0.1. Presumably, these common 
intersection points are just numerical coincidences, and will not occur in the general case 
(using other choices of e and D, for instance). However, the statistical effort for the data 
in Fig. [3] is rather substantial, and hence no such systematic parameter variation has been 
attempted. 

Figure [3b shows that the effective interfacial width w extracted from the fit to Eq. ([8]) 
increases from about w ^ 1.5 near the critical point of the thin film (the estimation of thin 
film critical points is discussed in the following sections) to about w ~ 2.4 for r^Jj = 0.8. 



87j. This complicated 



However, it is important to recall that the width w of the interface in the "soft mode" 
phase depends on both ?7p and the total film thickness D |83l. 
behavior results because the "intrinsic interfacial profile" 
5^ 



waves 



51 



52 



54 




891] is broadened by capillary 



551], but the long-wavelength part of the capillary wave spectrum is 



suppressed by the effective interface potential [38 



391 ] caused by the walls. For short range 



forces due to the wa 
width is 



83 



84 



85 



Is, as occurring here, the corresponding prediction for the mean square 



87] 



w 



wl[l + 



uj7r/A D 



+ const 



(9) 



2 + 00 Wo 

Here, Wq is the "intrinsic width", which should be related to the correlation length C,h 
along the coexistence curve in the critical region, wq = 2,^b, while the wetting parameter 



uj [38, 



39 



40 



90 



91 



92 



93| for Ising-like systems is cu ~ 0.8 and the (unknown) constant 



83 



M, 



851 ] can 



due to the short wavelength cutoff needed in the capillary wave spectrum 

be neglected near the critical point of the bulk. The intrinsic width should then vary with 

r/p as 

^0 = MvUvlcrit - l)~^ ^ ^ 0.63, (10) 



10 

with an amplitude factor wq which is presumably in the range 0.2 < wq < 0.5 (it is not 
accurately known since an unambiguous separation of intrinsic width and capillary wave 



broadening is hardly possible in interfacial profiles [50|, l87|]). Since for the chosen values of 
r]l we have D :^ wq ioi D = 10 and {uj7t/A)/(2 + uj) ^ 0.224, we expect that w ^ lA97y/w^ 
in our case, i.e. w in Fig. [3b should increase with an exponent 1^/2. Disregarding the results 
for ?7p = 0.8 and rj^ = 0.85, which are too close to r/pcrit ^^^ hence unreliable due to finite 
size effects, we find that the remaining data for L = 120 can be nicely fitted to a critical 
power law with the expected exponent z//2 = 0.315 (see insert of Fig. ^). 

Thus, it clearly would be of interest to obtain reliable data close to the bulk critical point, 
but then much larger systems would be required, and this would require very substantial 
computer resources, that are not available to us. But we emphasize the fact that no singular 
behavior can be observed when at fixed D we vary rj^ throughout the bulk critical region, 
passing the critical point. As an example. Fig. H] shows density profiles for the case D = 10, 
rjc = 0.195, L = 40 and three values of rj^ close to rj^^^j-n [Eq. ([T])]. One sees that profiles for 
r/p slightly above ?7p ^rit ^^^ slightly below it are hardly distinct from each other, all changes 
with respect to 1]^ are very gradual. 

A very interesting property is the correlation function of the colloidal particles in the 
interfacial region, zq — w<z<zq + w (see Fig. [5]). If we were to consider an unconfined 
interface, the capillary wave fluctuations would cause a power law decay of these fluctuations. 
Due to the confinement, the interface feels an effective potential, and this leads to the 
existence of a finite correlation length ^n of interfacial fluctuations, as discussed extensively 



in the literature 



m 



83 



84 



85 



871] . In simulations of a model for a symmetrical polymer 
mixture confined between competing walls, this correlation length was studied as a function 
of film thickness. Here we rather study this quantity as the interface localization transition 
is approached. Figure [5^ shows that the radial distribution function of colloidal particles in 
the interfacial regions is well described by the formula 

gdr) = const exp{—r/C,\\)/y/r . (11) 

Equation ( fTTI) was also shown to work very well in the case of the symmetric polymer 
mixture [831]. When rj^ approaches the value ?7p_crit(-^)5 o'^^ sees a strong increase of ^||, 
reflecting the expected critical divergence of ^|[ at the interface localization transition [which 
occurs at about ^pcrit(-^) ~ 1-13 ± 0.03]. Arguments have been given to show that for large 
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enough D there is a region of mean-field hke behavior, where ^n oc (1 — ''7p/''7pcrit(-^))~'^" 
with v\\ = 1/2, while very close toj?p^j,j^(D) the critical behavior should fall in the class of 
the two-dimensional Ising model 7^, u^ = 1. However, the accuracy of the data in Fig. [5b 
does not warrant an analysis of this crossover behavior. 



III. INTERFACE LOCALIZATION TRANSITION IN VERY THIN FILMS 



Following the procedures used in our earlier study of capillary condensation in the AO 
model, we carried out a finite size scaling analysis of the model with e = 3.0 for a slit pore 
which is only D = 3 colloid diameters thick. Varying the chemical potential and applying 



successive umbrella sampling [79[jthe probability distribution P{ric) is recorded. Applying 



suitable re- weighting techniques 9J] , one can apply the equal area rule 95|, l96|] to determine 
the chemical potential /^coex where the peak of -P(?7c) representing the vapor-like phase and 
the peak representing the liquid- like phase have equal weight. Figure Et shows typical data 
near the second order interface localization transition of the thin film, and Fig. [6b shows 
the fourth order cumulant f/4 as a function of t]^ for various L from L = 15 to L = 30. 
Introducing an order parameter M as M = rj^ — {i]c), the moments (M^) are defined as 

1 







(12) 



and U4 then is given as the ratio of the square of the second moment and the fourth moment. 



U^ = {My/{M^) 



(13) 



For large enough L, when finite size scaling [80|, [81|, |82] holds, a convenient recipe to find 
the critical point ^Jj ^rit is to record f/4 for different choices of L versus 1]^ tuning fi such that 
/^ = A^coexl^/p) and look for a common intersection point [80]. For r/p < r/pcrit one fixes /i by 
the criterion that (M^) is maximal {for ?7p > rjl, ^lit ^^is criterion is an alternative way to 
estimate /icoex(?7p)}- 

Figure [6^1 indicates the gradual change from a double peak distribution to a single peak 
distribution, which is a characteristic behavior for all second order phase transitions. Note 



that ?7p crit does not correspond to the value of r^p where P{ric) becomes fiat over a 



of rjc'. rather rj^^ ^^.^^^{D) still corresponds to a double peak distribution [80 



81 



Droad range 



82|. Figure 
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inb yields ^pcrit(-^ = 3) = 1.300 ± 0.005, i.e. a value very far away from r/p^j.;^. in the bulk 
[cf. Eq. ([1])]. Although it is somewhat disappointing that one cannot really find a unique 
intersection point of the cumulants f/4(?7p) for the various choices of L, one must recognize 
that for high enough resolution of the coordinate axes such a scatter is quite expected, due 
to residual corrections to finite size scaling 80(], and due to the statistical errors of the 



Monte Carlo data 



97l |. More disturbing is the fact that the cumulant intersections occur in 



a range of values in between the universal constants U* (2dim) and U* (3dim) for the two- 



and three-dimensional Ising model 98|, |99|], respectively. 



t/*(2dim) ^ 0.856, t/*(3dim) = 0.629 . 



(14) 



As Fig. [6b shows, intersections occur in the range 0.73 < U* < 0.80 (although there is some 



tendency of the intersection points to move upward with increasing L). On the other 



land, 



the slope of the cumulants at the intersection point, which is predicted to scale as SOj 

dUjd7]l oc L^/^ (15) 

yields an effective exponent rather close to the prediction u = 1 ioi the two-dimensional 
Ising model. 

Figure [7H shows simulation results for the order parameter m = [r]^ — vD/'^-i where the 
volume fractions of colloids r/f , rj^ are not read off from the peak positions in Fig. [6^, since 
for shallow peaks this would be a somewhat arbitrary procedure, but rather we take m as 
the first moment of the absolute value (|M|). Similarly, Fig. [7)d shows the "susceptibility" 
Xo = L^-D((M^) — (|M|)^). Both quantities are very strongly affected by finite size effects: 
Rather than exhibiting a power law decay, m oc (1 — ifp^critlVpY with (3 = 1/8, one finds that 
approaching r/p ^rit from above, the curves for m for the different values of L splay out and 



develop very pronounced "finite size tails" [80|, l95| for t]^ < rj^^^j-a. At r] = ripcrit '^^^ finds 



that the data are compatible with a power law decay (insert to Fig. [7^) 

m*(L)=m(L,r^^,,,iJocL-^/'^. 



(16) 



According to the two-dimensional Ising model, one would expect /?/z/ = 1/8. However, the 



straight line in the insert of Fig. [Th- rather indicates an effective exponent o 
0.20 ± 0.02. Likewise, the susceptibility maxima, which should scale as 



80 



81 



W/k 



82| 



cff 



xr^oc 



L^/" 



(17) 
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with the two-dimensional Ising value being 7/// = 1.75, rather suggest an effective exponent 
(7/'^)eff = 1-60 ± 0.03. Very roughly, these exponents are compatible with the hyperscaling 



relation [34] ^jv + 2/3/z/ = 2. Using the quoted effective exponents l/z/efr, (/5/'^)eflf and 
(7/'^)eff, one finds that on a scaling plot, where the variable t = |^p.crit/''?p — 1| is rescaled 
with L^/'^ and m or x are rescaled with L^^'^ or L"'^/'^, one finds reasonable data collapse 
(Fig. [8]). Such a partial success of a finite size scaling analysis, i.e. good data collapse is only 
found when effective exponents are used that deviate somewhat from the theoretical values, 
las already been seen for interface localization-delocalization transitions in the Ising model 
73I . l74l | and hence these problems are not a surprise in the present case. 

IV. OVERVIEW OF THE PHASE BEHAVIOR 

We now describe some of our results for other film thicknesses D. In principle, the same 
type of analysis was carried out for D = 5 and D = 7 as well, but it turned out that 
the distribution P(?7c) for rjc > ?7c,crit(-D) becomes increasingly asymmetric when D gets 
larger (Fig. [9]). Also the cumulant intersections get spread out over a rather large range 
of r/p (Fig. [TOj) . and these intersection points lie even in a range that is below the three- 
dimensional Ising value, Eq. (fT4l) . We interpret this finding as an indication that with D 
getting larger an increasing fraction of the critical region falls into the region of mean-field 



like behavior, as was theoretically predicted [74 1. 



Also for fixed D the accuracy, with which 1]^ critic) can be estimated, clearly deteriorates 
when e increases (Fig. [TTj) . Note that data ioi D = 5 and £ = 1.0 were already given in our 
preliminary communication [68|, the choice £ = 1.0 corresponds to a capillary condensation- 
type behavior, however. 

Figure fT2k shows estimates for the phase diagrams for the interface localization transition 
for e = 3 and three choices of D, while Fig. [T2b shows analogous data for D = 5 but varying 
e, and Fig. [13] shows a plot of rip,^^^^{D = 5) vs. e. One sees that miscibility is enhanced if 
either D decreases, or e increases, or both. 

Finally we turn to the variation of ^p^rit with e for the choice D = 5 (Fig. [T3l) . As 
found from a self-consistent field calculation for a symmetrical polymer mixture confined 
between competing walls [77| , the minimum of the curve rj^^^^^^ does not occur for the case 
of symmetric walls {e = 0), but for an asymmetric situation. It also is remarkable and 
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unexpected, that for large e the curve for rf^^^^i^ does not level off. 

Figure [H] shows the counterpart of the schematic Fig. [1] (left part), presenting in the 
plane of variables /ic and i]p{D) the numerical results for the coexistence curves between 
colloid-rich and polymer-rich phases, for the case of 2.0 < £ < 4.0, i.e. the region where 
interface-localization transitions occur (which are highlighted in the diagram by arrows). 
Note that unlike Fig. [H fic{oo) was not subtracted from /ic, thus the bulk coexistence is not 
simply the ordinate axis as in Fig. [1], but rather a nontrivial curve (which actually is not very 
different from a straight line). While for e = 2.0 there is still a small but systematic offset 
between the curves Hc{jf^,D = 5) and /ic(oo), for e = 3.0 and e = 4.0 the offset is almost 
negligibly small. The part of the curves /icoex(^p(-D)) to the left of ?7p^crit(-0) represents the 
state BI in the schematic phase diagram. Fig. [H where a delocalized interface occurs in 
the center of the film, separating the colloid-rich phase adjacent to the left wall and the 
polymer-rich phase adjacent to the right wall. 

At this point, we return to the density profiles at phase coexistence, and compare them 
for the same choice of rjl^ and D = 10, but different values oi e, e = 2.0 and e = 4.0 (Fig. [T5ll . 
For rjp = 1.5, the vapor-like phase reaches the same polymer density for both choices of e; 
the main difference concerns the colloid-rich side of the systems, the colloid enrichment at 
the hard wall is more pronounced for e = 4.0 than for e = 2.0. However, in the liquid- 
like colloid-rich phase the behavior is just the other way round: the layered profiles of the 
colloid-rich phase near the hard all are virtually identical, while the polymer enrichment 
near the right wall is more pronounced for e = 4.0 than for e = 2.0. When one studies the 
effect of varying e in the one phase region for r^p < r/p critl-^) however, one sees only a minor 
effect of e on the segregated structure where an interface has formed parallel to the walls 
(Fig. [T5b andfTSH). in particular for not extremely thin films. 

V. CONCLUSIONS 

In this paper, the Asakura-Oosawa (AO) model for colloid-polymer mixtures for a size 
ratio of polymers to colloids q = 0.8 was studied by Monte Carlo simulation, considering thin 
films of thickness D = 3 to D = 10 colloid diameters and confinement between asymmetric 
walls. One wall is simply a repulsive hard wall, to which the colloidal particles are attracted 
via depletion forces; the other wall exerts a square-well-type repulsive interaction (of the 
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range of the colloid ratio, and variable strength e = 0.5 to 4.0, in units of fc^T = 1.0). This 
study complements our earlier work on the AO model in the bulk, and under confinement 
between two equivalent hard walls, where capillary-condensation like phenomena occur; for 
the present model, we can smoothly interpolate from capillary condensation-like behavior 
for small e (e.g. e = 0.5 or 1.0), when both walls show some (though unequal) surface 
enrichment of colloids, to interface localization-type transitions, occurring for large e (e.g. for 
e varying from e = 2.5 to £ = 4.0). In the latter case, only the hard wall attracts colloids 
while the other wall attracts polymers. In this region, for large D the precise value of e 
has little effect on the observed density profiles. When one then increases the polymer 
reservoir packing fraction rj^ (which plays an analogous role as the inverse temperature does 
for thermally driven phase separation in small molecules mixtures), one observes that the 
enrichment layers of colloids and polymers at the walls gradually transform into two domains 
of coexisting colloid-rich and polymer-rich phases, separated by an interface parallel to the 
confining walls. We find that the temperature dependence of the width of this interface 
is considerably weaker than that of the bulk correlation length (or "intrinsic" interfacial 
width, respectively), and account for this finding in terms of capillary wave broadening of 
the interface. However, since for Z^ < 10 the interface profiles are strongly affected by 
layering of colloids near the hard wall, study of this broadening is difficult. 

Only far away from the bulk critical point can a sharp phase transition be observed, which 
we analyze by finite size scaling methods. While for D = 3 and not too large e the critical 
value ^p_crit(-^) can be rather accurately determined, and evidence can be found that the 
critical behavior falls in the universality class of the two-dimensional Ising model, for larger D 
and/or larger e the Monte Carlo data are strongly affected by problems of crossover between 
different universality classes and, thus, "/^pcritl-^) '^^^ t)e only estimated with rather modest 
accuracy, allowing no firm statements about critical exponents. Approaching the transition 
from r^p < rip^j.^^^{D), we find a strong increase of the correlation length ^|| describing the 
correlation of interfacial fiuctuations, but again the accuracy of our results would not suffice 
to estimate the value of the associated critical exponent. In view of the fact that even for 
the simple Ising model confined between competing boundaries a clarification of the critical 
behavior turned out to be very difficult, the problems encountered for the present more 
complicated model, which is strongly asymmetric even in the bulk, are not at all surprising. 

The fact that observation of interface localization does not require very special conditions 
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at the walls, but occurs for a broad parameter range, is encouraging for possible experimental 
tests of our results. We suggest that a repulsive interaction acting only on the colloids could 
be realized by creating a wall with a polymer brush at low grafting density. 

A very interesting problem, not accessible to the present grand-canonical Monte Carlo 
study, would be the dynamics of phase separation in such a confined thin film. We hope to 
report on such studies of a related model in the future. 

Acknowledgments: This work was supported in part by the Deutsche Forschungsgemein- 
schaft, SFB TR6/A5 and C3. 
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FIG. 1: Schematic phase diagrams of a colloid-polymer mixture confined between two parallel walls 
a distance D apart, in the grand-canonical ensemble where the polymer reservoir packing fraction 
rf is used as ordinate axis and the difference between the chemical potential of the colloids at 
bulk phase coexistence /icocx(-D = cxd) and the actual chemical potential of the colloids is used as 
abscissa (upper part). Phase coexistence in the bulk occurs along the vertical straight lines at 
Mcoex(c>o) — /x = 0. The left phase diagram refers to the case of symmetric walls, the right one 
to asymmetric walls. The lower part of the figure indicates the phases that occur in these phase 
diagrams (shaded regions denote colloid-rich domains): In the case of symmetric walls, a colloid- 
rich phase (AI) coexists along the line [i = ^cocx(-D, rf^) with a colloid-poor phase (All). In the case 
of asymmetric walls, the analogous phases are Blla, Bllb; Blla differs from Aj by the presence of 
a layer at the right walls where polymers are enriched, and Bllb differs from All by the fact that a 
surface enrichment layer of colloids exists only at the left wall. Finally, a state with a delocalized 
interface between colloid-rich and polymer-rich phases (BI) exists along the continuation of the 
line ^ = fJ-coex{D , rjp) beyond the critical point (^pcrit(^))' dotted line. The transition from BI to 
either Blla or Blla when one moves along the dotted line is termed interface localization transition. 
For further explanations cf. text. 
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FIG. 2: Colloid concentration profiles ?7c(-z) and polymer concentration profiles rjp^z) as a function 
of z for a thin film with asymmetric walls. At z = there is a hard wall, where the potentials 
U^^^{z) and t/w^p(-z) [Eqs. ([2|), ([3])] act on colloids (c) and polymers (p). Ki z = D there is another 
hard wall for both types of particles, with an additional square well repulsion acting on the colloids 
only [Eqs. (l3|), (H]), ([Sa]) . ([Fb]) . ([5c|) ] with a strength e = 2.5. Profiles were obtained at t/c = 0.18, 
ril = 0.70 (a), 7?c = 0.18, rf^ = 0.95 (b), r/c = 0.05, t?;^ = 1.20 (c), and r/c = 0.33, rf^ = 1.20 (d). For 
profiles (c) and (d), the choices t/c = 0.05 and 0.33 roughly correspond to the two branches of the 
coexistence curve in the bulk. 
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FIG. 3: Colloid density profiles ^c(-2) (a) and polymer density profiles ti^{z) (b) plotted vs. z for the 
case e = 3.0, D = 10, and choosing the average value of r/c equal to the diameter 5,5 = (ry^ + 7/^)/2. 
Here the origin at z = is put into the center of the slit pore. Broken curves denote fits with tanh 
profiles, as described in the text. Part (c) shows the dependence of the width w on rf . Critical 
values ff T-it in the bulk and in the thin film ?7pcrit(^) ^^'^ shown by arrows. Squares are for 
the choice D = 10, L = 120, circles for Z) = 10, L = 40. The insert shows a log-log plot of w 
vs. ?yp/?7p CTit ~ 1- The solid line in the insert is a fit with a power law with the exponent v/2 = 0.315 
(see text). 
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and 77c = 0.195. Three choices of r/p are included, as indicated. 
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FIG. 5: (a) Radial distribution function gc{r) of the colloidal particles, considering only distances 
r parallel to the walls, and particles confined in the interfacial region zq — w < z < zq + w, for 
the case L> = 10, L = 40, rjc = 0.195, e = 3 and two choices of ?yp, as indicated. Curves are fits 
to Eq. (jlip . (b) Plot of the parallel correlation length ^m extracted from fits as shown in part (a), 
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FIG. 6: (a) Logarithm lnP(?7c) plotted vs. t^c for the case -D = 3, L = 15, e = 3.0, and four values 
of r^p. (b) C/4 plotted vs. r^p, for the case D = 3, e = 3.0, and various choices of L, as indicated. The 
horizontal broken line indicates the range where intersections occur, 1.295 < ?/pcrit(^) — 1-305. 
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illustrates a slope l/i/n = 1.035. 
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FIG. 7: (a) Order parameter m plotted versus rf^ for D = 3, e = 3.0, and various choices of L 
as indicated in the figure. At rf ^^.j^ (indicated by a vertical dashed line) the resulting data are 
shown in a log-log plot versus L in the insert (the broken straight line indicates an effective exponent 
(/3/7)cflf ~ 0.2). (b) Susceptibility Xc plotted versus rf^ , for the same choices as in (a). Insert shows 
the maxima on a log-log plot vs. L (broken straight line indicates an exponent (7/j^)eflf ~ 1-6). 
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FIG. 8: Scaling plot of the order parameter (a) and the susceptibility (b), using the three largest 
values of L in Fig. [71 and the effective exponents (omitting here the index "eff " ) l/v = 1.035, iS/u = 
0.2, and j/u = 1.675. The straight lines indicate the theoretical slope of the scaling functions for 
large values of tV-'^ , namely (3 = 1/8 (a) and 7 = 7/4 (b). 
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FIG. 9: (a) Logarithm of the probability distribution of the colloid volume fraction, lnP(r/c), 
plotted vs. r/c for D = 5, L = 20, e = 3.0, and four choices of rj^, as indicated, (b) Same as (a), 
but for D = 7. 
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FIG. 10: Cumulant U^ plotted vs. rf^ for various choices of L as indicated, for e = 3.0 and the 
film thickness D = 5 (a) and D = 7 (b). Inserts show plots of InYi vs. InL, where Yi = dU^/drf . 
The broken straight lines indicate the effective exponent (l/i/)cff, with (a) {l/v)cs ~ 1.08 and (b) 
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FIG. 12: Phase diagrams in the plane of variables rf^ and r^c for the case L = 30, e = 3.0, and three 
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FIG. 13: Plot of r?;;^erit(^ = 5) vs. e. 
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FIG. 14: Plot of the colloid chemical potential versus rf . The coexistence curve for asymmetric 
thin films, for D = 5 and three values of e, e = 2.0,3.0 and 4.0, as indicated, are compared to the 
coexistence curve of the colloid-polymer mixture in the bulk (solid line). Arrows show the critical 
points ryp^j.j^(L'). Note that the bulk critical point (^pj.j.jt = 0.766) is far beyond the scale of the 
diagram. For ry!, < ripcnti-^) ^^^ meaning of the shown curves denotes the coexistence of domains 
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b) liquid, D=5,r|^=1.5 
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c) one-phase region, D = 5, r\[^= 1 .0 d) one-phase region, D= 10, r|'^=0.9 
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FIG. 15: Comparison of the density profiles of colloids and polymers in the coexisting vapor-like 
phase (a) and liquid-like phase (b) between the cases e = 2.0 and 4.0, for rf = 1.5 and D = 5. 
Panels c) and d) show the comparison of the density profiles of colloids and polymers between the 
cases e = 2.0 and 4.0, but for the one-phase region in the "soft mode" phase. Panel c) displays the 
profiles for D = 5 and r^Jj = 1.0, panel d) the ones for D = 10 and ?7p = 0.9. 
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